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In this work a symmetry of universal finite-size scaling functions under a certain anisotropic scale 
transformation is postulated. This transformation connects the properties of a finite two-dimensional 
system at criticality with generalized aspect ratio p > 1 to a system with p < 1. The symmetry 
is formulated within a finite-size scaling theory, and expressions for several universal amplitude 
ratios are derived. The predictions are confirmed within the exactly solvable weakly anisotropic 
two-dimensional Ising model and are checked within the two-dimensional dipolar in-plane Ising 
model using Monte Carlo simulations. This model shows a strongly anisotropic phase transition 
with different correlation length exponents v\\ / u± parallel and perpendicular to the spin axis. 



The theory of universal finite-size scaling (UFSS) func- 
tions is a key concept in modern understanding of con- 
tinuous phase transitions [EL 0, 01. In particular, it is 
known that the UFSS functions of a rectangular two- 
dimensional (2D) system of size L\\ x L± depend on the 
aspect ratio L\\/L± [f|. For instance, in isotropic systems 
the scaling function at criticality U c of the Binder cumu- 
lant U = 1 — i(rn 4 )/(m 2 ) 2 where (m n ) is the n-th 
moment of the order parameter, is known to be a univer- 
sal function U C (L\\/L±) for a given boundary condition. 
This quantity has been investigated by several authors 
in the isotropic 2D Ising model with periodic boundary 
conditions |^, Jj|, while the influence of other boundary 
conditions on U C (L\\/L±) has recently been studied in 
Refs. § §. 

In weakly anisotropic systems, where the couplings are 
anisotropic (Jn ^ J± in the 2D Ising case), the correla- 
tion length of the infinite system in direction p = II, ± 
becomes anisotropic and scales like ^°°' ) (t) ~ £/ t t _1/ near 
criticality. (t = (T — T c )/T c is the reduced temperature 
and we assume t > without loss of generality.) This 
leads to a correlation length amplitude ratio £n/£_i_ differ- 
ent from unity. The UFSS functions then depend on this 
ratio, i.e. U c = U c (Li\/L±,£\i/£±). However, isotropy can 
be restored asymptotically by an anisotropic scale trans- 
formation, where all lengths are rescaled with the cor- 
responding correlation length amplitudes £ M [ ^o[ O , fl2"|| . 
Thus the UFSS functions depend on L\\/L±_ and £n/£j_ 
only through the reduced aspect ratio (L\\/£\\)/(L±/£x)- 

In strongly anisotropic systems both the amplitudes £ M 
as well as the correlation length exponents are differ- 
ent and the correlation length in direction p scales like 



(i) 



Examples for strongly anisotropic phase transitions are 
Lifshitz points as present in the anisotropic next 
nearest neighbor Ising (ANNNI) model (Til Ej| pf, or 



the non-equilibrium phase transition in the driven lat- 
tice gas model [Ej], Furthermore, in dynamical sys- 
tems one can identify the ll-direction with time and the 
-L-direction(s) with space [Ej}), which in most cases give 



strongly anisotropic behavior. 

Using the same arguments as above we conclude that 
UFSS functions of strongly anisotropic systems depend 
on the generalized reduced aspect ratio (cf. fq[) 



p = L||L ± e /r 5 , with r$ = |n£] 



(2) 



being the generalized correlation length amplitude ratio, 
and with the anisotropy exponent 9 = v\\/v±_ Jl9). Up 
to now no attempts have been made to describe the de- 
pendency of UFSS functions like U c (p) on the shape p 
of strongly anisotropic systems. In particular, it is not 
known if the anisotropy exponent 9 can be determined 
from U c {p). This problem is addressed in this work. 

Consider a 2D strongly anisotropic finite system with 
periodic boundary conditions. When the critical point 
of the infinite system is approached from temperatures 
t > 0, the correlation lengths £ M in the different directions 
/i are limited by the direction in which from Eq. ([j]) 
reaches the system boundary first Qj. For a given volume 
N = L\\L± we define an "optimal" shape Ljj pt x L' 



which both correlation lengths £p reaches the system 
boundary simultaneously, i. e. 
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for some temperature t > (Fig. 0a). We immediately 
find using Eqs. (Q, ^|) that the optimal shape obeys p pt = 
1 for all N, giving Ljj pt = r^(L^ pt ) e . A system of optimal 
shape should show the strongest critical fluctuations for a 
given volume N as the critical correlation volume £ii, c £_l,c 
spans the whole system. 

At the optimal aspect ratio p = 1 the correlations are 
limited by both directions II and _L (Fig. [j]a). If the 
system is enlarged by a factor b > 1 in the ll-direction 
(Fig. 0b) , the correlation volume may relax into this di- 
rection but does not fill the whole system due to the 
limitation in ^-direction. A similar situation with ex- 
changed roles occurs if the system is enlarged by a fac- 
tor b > 1 in the _L-direction (Fig. [j]c). We now assume 
that systems (b) and (c) are similar in the scaling re- 
gion L° pt — > oo, i. e. that their correlation volumes are 
asymptotically equal. 
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Figure 1: Three systems with different aspect ratio p (Eq. ^f) 
at criticality. In (a) the critical correlation volume fu, c £_L,c 
(shaded area) spans the whole system, while in (b) and (c) 
correlations are limited by symmetric finite-size effects. 



Hence we can formulate a symmetry hypothesis: Con- 
sider a system with periodic boundary conditions and 
optimal aspect ratio p = 1 at the critical point. If this 
system is enlarged by a factor b > 1 in ll-direction, it be- 
haves asymptotically the same as if enlarged by the same 
factor b in _L-direction. 

To formulate this hypothesis within a finite-size scaling 
theory, we consider a 2D strongly anisotropic system of 
size L\\ x L± which fulfills the generalized hyperscaling 
relation 2 — a = i>\\+i'± ||. For our purpose it is sufficient 
to focus on the critical point. The universal finite-size 
scaling ansatz |], |j| |], |(| for the singular part of the 
free energy density f c = F S)C / (NksTc) reads [M 



f c (L ii ,L 1 _)r, b -^Y c (b ll ,b 1 _) 



(4) 



with the scaling variables 6 M = A^L^/^, where A is a 
free scaling parameter. The scaling function Y c is uni- 
versal for a given boundary condition, all non-universal 
properties are contained in the metric factors £ M . These 
metric factors occur due to the usual requirement that 
the relevant lengths are L^/Q, 00 \t) near criticality, and 
cannot be absorbed into A in contrast to isotropic sys- 
tems. For the three systems in Fig. |l| we set A = 
{LT/^r 1 '^ to get 
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(5b) 
(5c) 



The proposed symmetry hypothesis states that for b > 1 
Eqs. (5b) and (pq) are asymptotically equal in the scaling 
region where L° pt is large, 



Opts. 
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Hence the scaling function Y c has the simple symmetry 
^(6,1)^^(1,6). (7) 

To rewrite Y c as function of the generalized aspect ratio 
p (Eq. (||)) instead of the quantities 6^, we set b± = 1 in 

system (c) and get, as then A = {bL°f /£x)~ 1 ' vx , 



f c {L°?\bLT) 



OptN 



N 



Y c (b- e ,l) 



(8) 



Eqs. (|k]) and (||) are identical and we conclude that 
&y c (l, b) = b~ e Y c (b~ e , 1). At this point it is convenient to 
define the scaling function Y c (b) = bY c (b, 1) which fulfills 



/c(£ii,£l)~-^c(p). 



(9) 



For this scaling function the symmetry reads 

Y c (p) >2 Y c {p~ e ). (10) 

We see from Eq. @ that the critical free energy density 
f c is a universal function of the reduced aspect ratio p — 
LwLj 6 /r^ without any non-universal prefactor, and that 
at criticality all system specific properties are contained 
in the non-universal ratio from Eq. (||). 

Ansatz Eq. (|j) can also be made for the inverse spin- 
spin correlation length at criticality |2(j 



£-i(L H ,L ± )~-±Xt l , c (h,b ± ) 



(11) 



The proposed symmetry gives X^^b, 1) = Xp iC (X,b), 
where p denotes the direction perpendicular to p. Hence 
the correlation volumes £ll,c£-L,c of system (b) and (c) in 
Fig. [| are indeed equal as assumed above and become 
6i,ca j c~f^ c 1 (6,l)Xl i 1 c (6,l). 

The correlation length amplitudes A^ in cylindrical ge- 
ometry (6^ 



oo, bp, 



1), which can be calculated ex- 
actly for many isotropic two-dimensional models within 
the theory of conformal invariance |^l| generalize to the 
strongly anisotropic form H 

A*= lim Lr«W"A l im ^ C (L U ,L ± ). (12) 



Inserting Eq. (|11|) they become 



4 = ^cV,i), 4 



-1/6 



Xjia.oo) (13) 



which shows that in general A% ist not universal. The 
symmetry hypothesis states that both limits of the scal- 
ing function X MjC are equal. Denoting this universal limit 
A e := ^^(oo, 1) = X±* c (l,co) we obtain A^ = r ( A^ 

A^ as well as the amplitude relations 



and Aj- — r, l '' 



A\ +e 
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(14) 
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(oo) (i) = (logcoth(/3J M ) - 2/3J;,)- 1 
The amplitude ratio at the 



These predictions can be checked within the exactly 
solved weakly anisotropic 2D Ising model with different 
couplings Jn and J±, where the paramagnetic correla- 
tion length reads £ 
with (3 = l/k B T § 
critical point sinh(2/3 c Jn) sinh(2/j c J±) = 1 ||22j becomes 
rc — sinh(2/3 c Jn) (g3|. On the other hand, the inverse cor- 
relation length amplitudes in cylinder geometry Eq. ( |l^ ) 
has been calculated (24|] to give A£ = ^ sinh(2/3 c J p ), 
which immediately yields Eqs. ([l3]) if we insert the well 
known universal value = 4/ir p5| . The left re- 
lation of Eqs. ( |l4| ) has already been derived for sev- 
eral weakly anisotropic models, where it simplifies to 
i^Af @, Eq. (7)]. 
To check the symmetry numerically in strongly an- 
isotropic systems, we now focus on the Binder cumulant 
U. The scaling ansatz at criticality Eq. (0) becomes 
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U c (L n ,L ± ) 
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■U c (h,b x ) = tJ c (p) 



(15) 



with the scaling function U c (b) = U c (b, l)/b, and the cal- 
culation is completely analogous to the free energy case. 
The symmetry hypothesis for the cumulant scaling func- 
tions U c and U c thus reads (cf. Eqs. 



(16) 



u c (b,i) b ^u c (i,b), u c ( P ) p =u c { P - e ). 

The generalization of the cumulant amplitude A 1 /, 
to strongly anisotropic systems is similar to Eq. (TL2 
gives 



A u = T lim L n 

Lu — >oc 



lim Z, M [/ c (I,||,LjJ. 



(17) 



Inserting the scaling ansatz Eq. ( |15[ ) we now find 

^=r e L> c (oo,l), 



A} 



C/ c (l,oo), (18) 



which again are in general not universal. Using the 
symmetry hypothesis we can define Ajj := U c {oo, 1) 
f7 c (l,oo) and get A v ^nAu, A\ 
the identities (cf. Eqs. (Q) 



' Au as well as 



A v +e = A(Au) 6 



,4-1 



(19) 



The cumulant scaling function U c (p) must be extremal 
at p = 1 due to symmetry. Furthermore, as a deviation 
from the optimal aspect ratio p — 1 reduces the cumu- 
lant, it has a maximum at this point A sketch of 
U c {p) for an assumed anisotropy exponent 9 — 2 is de- 
picted in Fig. |. For p > 1 both C4(p) and U c (p' = p^ ) 
collapse onto a single curve, reflecting the proposed sym- 
metry. It is obvious from Fig. ^ that U c (p) (and thus 
also Y c (p) from Eq. (|i"o|)) can not be analytic at p = 1 in 
strongly anisotropic systems, as the two branches U c (p) 
and Uc(p') identical for p > 1 fork at p = 1 [M. On 
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Figure 2: Sketch of critical cumulant scaling functions U c (p) 
and Uc(p') with p = p~° for assumed anisotropy exponent 
= 2. We have U c {p>l) ~ Avjp and U c {p<X) ~ Ayp 179 . 
For p > 1 C7 c (p) fulfills (7 C (p) = £7 c (p'). 



the other hand, Y c {p) and t7 c (p) can be analytic at p = 1 
if the anisotropy exponent 6 = 1, as in the case of the 
isotropic 2D Ising model Eq. 3.37]. 

To check the symmetry hypothesis in a strongly 
anisotropic system, I performed Monte Carlo simulations 
of the two-dimensional dipolar in-plane Ising model pcj 



(ij) *7%7 



2(^) 2 



(20) 



with spin variables a = ±1, ferromagnetic nearest neigh- 
bor exchange interaction J > 0, and dipole interaction 



lu > 0. The distance ry = (? 



) between spin cr^ and 



cTj is decomposed into contributions parallel and perpen- 
dicular to the spin axis. In the simulations the Wolff 
cluster algorithm |2q| for long range systems proposed by 
Luijten and Blote ||29| was used, modified to anisotropic 
interactions. In contrast to earlier work [^o[ ||l[] using 
renormalization group technics it is found that this model 
shows a strongly anisotropic phase transition. The de- 
tails of the simulations will be published elsewhere |^o|| . 

After T c was determined, systems with constant vol- 
ume N = L\\L± were simulated, which was chosen to 
have a large number of divisors in order to get many dif- 
ferent aspect ratios (e. g. N = 2 6 3 3 5 2 = 43200 has 84 
divisors). The resulting critical cumulant U C (L\\L^ ) for 
two different volumes N — 4320, 43200 is depicted in the 
inset of Fig. [|. As expected, both curves have the same 
maximum value U c (l) — 0.555(5) at criticality. With 
variation of 9 the curves are shifted horizontally and col- 
lapse for 9 = 2.1(3), with maximum at — 0.415(40). 
To check the proposed symmetry we fold the left branch 
with p < 1 (open symbols) to the right and rescale the 
p-axis with 9. The resulting data collapse for p > 1 is 
shown in Fig. @. This collapse and the additional condi- 
tion that both curves must go to zero as Au/p allows a 
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Figure 3: Cumulant U C (L\\, L±_) of the dipolar in-plane Ising 
model (Eq. (po|)) for dipole strength to/ J = 0.1 and system 
size N = 43200 at the critical point k B T c /J = 2.764(1). The 
data points collapse for p > 1 if we set 9 = 2.1(3) and = 
0.415(40), giving the universal amplitudes t7 c ( 1) = 0.555(5) 
and Ay = 3.5(2). The inset shows U c as function of the un- 
reduced generalized aspect ratio L\\LJ° for system size N = 
43200 (circles) and N = 4320 (triangles). 



precise determination of 6 and as well as of the uni- 
versal amplitude Ajj = 3.5(2). 

In conclusion, I postulate a symmetry of universal 
finite-size scaling functions under a certain anisotropic 
scale transformation and generalize the Privman-Fisher 
equations Jl| to strongly anisotropic phase transitions on 
rectangular lattices at criticality. It turns out that for 
a given boundary condition the only relevant variable is 
the generalized reduced aspect ratio p = L\\L e ± /r£ and 
that e. g. the free energy scaling function Eq. @ obeys 
the symmetry Y c (p) "= Y c (p~ e ). At criticality, the free 
energy density f c , the inverse correlation lengths £ M , C , 
and the Binder cumulant U c are universal functions of 
p, without a non-universal prefactor. All system specific 
properties are contained in the non-universal correlation 
length amplitude ratio (Eq. Q)). 

The generalization to higher dimensions is straightfor- 
ward |^(|, an interesting application would be the precise 
determination of the exponent 9 at the Lifshitz point 
of the three-dimensional ANNNI model [§]| An 
open question is the validity of the proposed symmetry 
in non-equilibrium systems with appropriate boundary 
conditions, which recently have been shown to exhibit 
Privman-Fisher universality ||. 
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